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Abstract 

Introduction: There is a major need to better understand the molecular basis of triple negative breast cancer 
(TNBC) in order to develop effective therapeutic strategies. Using gene expression data from 587 TNBC patients we 
previously identified six subtypes of the disease, among which a mesenchymal-stem like (MSL) subtype. The MSL 
subtype has significantly higher expression of the transforming growth factor beta (TGF-(3) pathway-associated 
genes relative to other subtypes, including the TGF-(3 receptor type III (T(3RIII). We hypothesize that T(3RIII is tumor 
promoter in mesenchymal-stem like TNBC cells. 

Methods: Representative MSL cell lines SUM159, MDA-MB-231 and MDA-MB-157 were used to study the roles of 
T(3RIII in the MSL subtype. We stably expressed short hairpin RNAs specific to T(3RI 1 1 (T(3RIII-KD). These cells were then 
used for xenograft tumor studies in vivo; and migration, invasion, proliferation and three dimensional culture studies 
in vitro. Furthermore, we utilized human gene expression datasets to examine T(3RIII expression patterns across all 
TNBC subtypes. 

Results: T(3RI 1 1 was the most differentially expressed TGF-(3 signaling gene in the MSL subtype. Silencing T(3 Rl 1 1 
expression in MSL cell lines significantly decreased cell motility and invasion. In addition, when T(3RIII-KD cells were 
grown in a three dimensional (3D) culture system or nude mice, there was a loss of invasive protrusions and a 
significant decrease in xenograft tumor growth, respectively. In pursuit of the mechanistic underpinnings for the 
observed T|3RIII-dependent phenotypes, we discovered that integrin-a2 was expressed at higher level in MSL cells 
after T|3RIII-KD. Stable knockdown of integrin-a2 in T|3RIII-KD MSL cells rescued the ability of the MSL cells to migrate 
and invade at the same level as MSL control cells. 

Conclusions: We have found that T|3RIII is required for migration and invasion in vitro and xenograft growth in vivo. 
We also show that T|3RIII-KD elevates expression of integrin-a2, which is required for the reduced migration and 
invasion, as determined by siRNA knockdown studies of both T|3RIII and integrin-a2. Overall, our results indicate a 
potential mechanism in which T|3RIII modulates integrin-a2 expression to effect MSL cell migration, invasion, and 
tumorigenicity. 
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Introduction 

The term 'triple negative breast cancer' (TNBC) is used to 
classify the 10% to 20% of all breast cancers that lack estro- 
gen receptor (ER) and progesterone receptor (PR) expres- 
sion as well as amplification of the human epidermal 
growth factor receptor 2 (HER2) [1]. Disease heterogeneity 
and the absence of well-defined molecular targets have 
made treatment of TNBC challenging. There is a major 
need to understand better the molecular basis of this type 
of breast cancer in order to develop effective therapeutic 
strategies [1]. In a previous study, we performed gene 
expression (GE) analyses and identified six distinct mo- 
lecular TNBC subtypes with unique biological drivers 
[2], including one that was enriched for mesenchymal- 
associated genes termed mesenchymal-stem like (MSL). 
The MSL subtype is characterized by increased expres- 
sion of genes related to transforming growth factor beta 
(TGF-p) signaling as well as pathways that play roles in 
extracellular matrix (ECM), focal adhesion, cell motility 
and cell differentiation [2]. Of note, TGF-p receptor 
type III (TpRIII) (gene symbol: TGFBR3) was among the 
differentially expressed TGF-p signaling components in 
the MSL subtype. 

The TGF-p signaling pathway has been implicated in 
cancer initiation and progression through tumor cell au- 
tonomous and non-autonomous signaling [3,4]. Initially 
identified as a tumor suppressor and then as a mediator 
of tumor progression, TGF-p signaling demonstrates di- 
verse capabilities in cancer. The TGF-p pathway suppresses 
tumor growth through regulation of epithelial and stromal 
cell signaling [5]. Dysfunction of the pathway leads to car- 
cinoma progression and metastasis [3]. While there has 
been significant focus on TGF-p receptor type I (TpRI) and 
TGF-p receptor type II (TpRII), research on TpRIII has 
lagged. Prior studies have demonstrated that TpRIII can 
regulate TGF-p signaling either via delivering TGF-p2 
ligand to TpRII [6-9] or by binding to the cytoplasmic 
domain of TpRII, forming an active TpRI-TpRII signal- 
ing complex [10-13]. Currently, analysis of gene expres- 
sion data sets generated from multiple cancer types 
indicates that TpRIII expression is decreased in higher- 
grade cancers [14-17]. However, the role of TpRIII is 
controversial in breast cancer, since it has been reported 
that TpRIII can act as either a tumor suppressor or pro- 
moter in this cancer [18,19]. 

In the current study, we focused our investigations on 
the functional role of TpRIII in the MSL subtype of 
TNBC. We used a loss-of-function approach in represen- 
tative MSL cell lines to demonstrate that TpRIII is re- 
quired for maintenance of tumorigenicity in MSL TNBC 
cell lines and that regulation of integrin-a2 (gene symbol: 
ITGA2) expression is mechanistically involved in the ob- 
served pheno types. This study demonstrates that TpRIII 
promotes the in vivo growth of a subset of TNBC and 



provides a pre-clinical rationale for consideration of TpRIII 
as a potential target for further discovery efforts. 

Materials and Methods 

Cell culture 

SUM159 cells (Asterand, Detroit, MI, USA) were main- 
tained in (Dulbeccos) Modified Eagle s Medium: Nutri- 
ent Mixture F12 ((D)MEM-F12, GIBCO, Grand Island, 
NY, USA) supplemented with 5% fetal bovine serum 
(FBS) (GIBCO) and 0.5 ug/ml hydrocortisone. MDA- 
MB-231 and MDA-MB-157 (ATTC, Manassas, VA, 
USA) were maintained in (D)MEM (GIBCO) supple- 
mented with 10% FBS. Stable TpRIII-KD SUM159 cell 
lines were generated by lentiviral infection with virus 
carrying four independent short hairpin RNA (shRNA) 
clones (sequence-verified shRNA, pLKO.l-puro), 
(Sigma- Aldrich, St. Louis, MO, USA), Mission shRNA li- 
brary #SHCLNG-NM_003243: clone#TRCN0000033433 
(TpRIII-KD), clone#TRCN0000359000 (TpRIII-KD2), 
clone#TRCN0000359001 (TpRIII-KD3), and clone# 
TRCN0000359081 (TpRIII-KD4)) followed bypuromycin 
selection (Invitrogen-Life Technology, Inc, Carlsbad, CA, 
USA). MDA-MB-231 and MDA-MB-157 were stably in- 
fected with clone# TRCN0000033433. Integrin-a2 was 
stably knocked down in TpRIII-KD MSL cells using 
lentiviral particles carrying shRNA to integrin-a2 (a2- 
KD) (Sigma-Aldrich, Mission shRNA validated library, 
#SHCLNG-NM_002203, clone#TRCN0000308081). 

Three-dimensional culture assay 

The wells in 48-well plates were coated with 50 \A of growth 
factor reduced BD Matrigel (BD Biosciences #356231, 
San Jose, CA, USA) and allowed to polymerize at 37°C 
for 15 minutes. Then, 5 x 10 5 cells were resuspended in 
200 \A of growth factor reduced BD Matrigel and plated 
onto the matigel- coated wells. Plates were incubated for 
30 minutes after which 1 ml of media was added to the 
top of the matrigel. Media was replenished every 48 hours. 
Images were taken at day six. Quantification of the images 
was performed using Fiji Software. 

Cell proliferation assays 
Cell counts 

Cells were plated into six-well plates at a density of 1.25 x 
10 5 cells/ well. The following day cells were treated with 
1 ng/ml TGF-pl (R&D Systems, #102-B1, Minneapolis, 
MN, USA) and TGF-P2 (R&D Systems, #102-B2). After 
72 hours treatment with TGF-p, viable cells were counted. 

3 H-Thymidine incorporation assay 

A total of 2.5 x 10 4 cells were plated in a 24-well dish 
and allowed to grow overnight. The next day the medium 
was aspirated and replaced with complete medium con- 
taining +/-TGF-pl or TGF-p2 (1 ng/ml). The cells were 
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then subjected to [ 3 H] thymidine incorporation assay as 
previously described [20] . 

Migration and invasion assays 

Magnetic attachable stencils migration assays 

This migration method serves as a more reproducible alter- 
native to the scratch assay. The use of magnetic force to at- 
tach stencil to the multi-well plates is a new strategy that 
creates defined and reproducible cell-free voids for quanti- 
tation of cell migration and has been well characterized and 
described by Ashby et al [21]. Magnetic attachable stencils 
(MAtS) were attached to the surfaces of each well of a 12- 
well plate by placing a platform with magnets underneath 
and in direct contact with the 12- well plate. Cells were then 
plated in triplicate at 7.5 x 10 5 cells per well around the 
MAtS in serum-free media. The next day the MAtS were 
removed and cells were treated with 1 ng/ml TGF-pl 
(R&D Systems, #102-B1) and 1 ng/ml TGF-(32 (R&D Sys- 
tems, #102-B2). Gap closure was quantified (Tscratch soft- 
ware) at both 0 and 24 hours and percent of closure 
determined with the following equation: percent of closure 
= average of ((gap area: 0 hour) - (gap area: 24 hours))/(gap 
area: 0 hour) using images from 12 different microscopic 
fields per well (4X magnification). 

Transwell assays 

Migrations (Costar, #3422, Tewksbury, MA, USA) were 
conducted by plating 2.5 x 10 4 cells in the top of the trans- 
well and media with 10% FBS in the bottom of the well to 
act as a chemoattractant. Cells were fixed in 4% paraformal- 
dehyde and stained with 4', 6-diamidino-2-phenylindole 
(DAPI). Quantification was performed by taking pictures of 
multiple regions of the membrane after which cells' nuclei 
were counted using Metamorph software. The same migra- 
tion assay was used to measure blocked integrin-a2 func- 
tion. The TpRIII-KD cells were incubated for 30 minutes 
with integrin-a2 blocking antibody (Abeam, #ab24697, 
Cambridge, MA, USA) washed two times with PBS and 
plated in the top of the transwell. Invasion assays were 
conducted by plating 5 x 10 5 cells in serum-free media 
in the upper chamber, pre-coated with growth factor re- 
duced matrigel. In the bottom chamber (D)MEM with 
10% FBS was used as a chemoattractant (BD Biosci- 
ences, #354483). Cells that had invaded through the 
matrigel were fixed in 4% paraformaldehyde and stained 
using DAPI. Quantification of cells that invaded into the 
matrigel was performed using the same protocol as de- 
scribed for the transwell assays. 

Xenograft tumor studies 

One milllion cells embeded in collagen were implanted 
into the number four gland of six- to eight-week-old 
female athymic nude- Foxnl nu/nu mice (purchased from 
Harlan Sprague- Dawley, Inc., Indianapolis, IN, USA). 



Mice were monitored weekly for tumor growth. Tumor 
measurements were performed once a week for five 
weeks after palpable tumors formed. Tumor volume was 
measured at the indicated times with calipers, and tumor 
volumes were calculated as width 2 x length/2. All mouse 
experiments were approved by the Vanderbilt University 
Institutional Animal Care and Use Committee (IACUC). 

Luciferase reporter assay 

Cells were seeded at a density of 2 X 10 4 cells/well in 12- 
well tissue culture plates. The following day, the cells were 
transiently transfected using Transfectin lipid reagent fol- 
lowing the manufacturers protocol (Bio-Rad #170-3351, 
Hercules, CA, USA). Cells were transfected with 1.5 \ig 
3TP-Lux [22] or CAGA(9)-Luc [23]. pRL-CMV-renilla 
(Promega #E2261, Madison, WI, USA) was co-transfected 
and used as an internal control to correct for transfection 
efficiency. Eighteen hours after transfection, cells were 
treated with 1 ng/ml TGF-(3l or TGFp-2 (R&D Systems, 
#102-B1 and #102-B2, respectively). Twenty-four hours 
after TGF-p treatment, cells were harvested and assayed for 
promoter specific luciferase activity using a Dual-Luciferase 
Reporter Assay System (Promega #E1910) according to the 
manufacturer s protocol. Luciferase activity was measured 
using a BD/Pharmigen Monolight 3010 luminometer. 

RNA preparation and quantitative PCR 

RNA was isolated and purified using an RNeasy Mini 
Kit and an RNase-Free DNase Set (Qiagen, Valencia, CA, 
USA). A total of 750 \i% of RNA was used to synthesize 
cDNA using Superscript III reverse transcriptase as de- 
scribed by the manufacturer (Invitrogen). Bio-Rad iCy- 
cler and CFX96 machines were used for qPCR employing 
Power SYBR Green (Applied Biosystems, Carlsbad, CA, 
USA) or SsoAdvanced SYBR Green Supermix (Bio-Rad), 
respectively. C t values were normalized to GAPDH for 
statistical analyses. Primer sequences are available in 
Additional file 1. 

Immunoblotting 

Standard protein preparation and electrophoresis proce- 
dures were used as described [4]. Western membranes 
were blocked in 5% milk and incubated with primary anti- 
body overnight. The antibody list with concentrations and 
the catalog numbers are available in Additional file 1. 

Flow cytometry 

Cells were detached using Accutase (Life Technologies), 
pelleted, washed and counted. One million cells were incu- 
bated with TpRIII antibody (Cell Signaling, #5544, Danvers, 
MA, USA) for 30 minutes, washed, and then incubated at 
4°C with Alexa Fluor 488 conjugated secondary antibody 
(Life Techologies, #A11034) for 30 minutes. One million 
cells were labeled with fluorescence-conjugated integrin-a2 
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antibody (BioLegend, #314308, San Diego, CA, USA) for 
30 minutes at 4°C. Cells were washed three times then ana- 
lyzed on a FACSCalibur flow cytometer (Becton Dickinson, 
San Jose, CA, USA) using CellQuest Pro software. Data 
were analyzed with Flowjo software (Tree Star). 

Microarray gene expression analysis 
Public database analysis 

Human tissue and cell line microarray datasets were an- 
alyzed using GeneSpring GX 12.0 microarray analysis 
software (Agilent). Previously published TNBC gene ex- 
pression profiles (n = 587 patients) [2] consisting of pub- 
licly available microarray data sets (the GEO registration 
numbers are referenced in Additional file 1) were obtained 
and processed as previously described and were in compli- 
ance with ethical requirements [2]. Comparisons between 
expression of TGFBR3 and ITGA2 for different TNBC 
subtypes were performed in R 3.0.1 [24] using the t test 
function for paired two-tailed Students £-tests and graph- 
ically represented using ggplot2 [25]. 

In vitro three-dimensional culture analysis 

vRNA was extracted from SUM159 three-dimensional cul- 
ture samples and hybridized to the human gene LOST 
array, scanned with Affymetrix using AGCC v. 3.2.4 and 
then analyzed in R 3.0.1 using the oligo package. Samples 
were normalized with the RMA algorithm, genes were an- 
notated with the pd.hugene.l.O.stvl package, and differen- 
tial gene expression analysis was conducted using the 
limma package. The three-dimensional culture microarray 
data discussed in this publication have been deposited in 
the National Center for Biotechnology Information 
(NCBI) s Gene Expression Omnibus [26] and are accessible 
through GEO Series accession number GSE54756 [27]. 

Statistical analysis 

All data were analyzed using the unpaired two-tailed 
Students t test (GraphPad Prism 5). Error bars show 
mean ± SEM. A two-sided P value less than 0.05 was 
considered significantly different. 

Results 

Human mesenchymal stem-like triple negative breast 
tumors and representative cell lines have increased 
TPRIII expression 

Using a gene expression data set generated from 587 TNBC 
tumors, we examined the relative mRNA levels of TGF-p 
receptors and ligands across subtypes of TNBC. We ob- 
served elevated expression of TGFBR3 in basal-likel (BL1), 
mesenchymal (M) and MSL tumors (Figure 1A). The high- 
est relative level of TGFBR3 expression was in the MSL 
subtype (Figure IB). Average probe intensities for the TGF- 
p receptors I and II as well as TGF p ligands 1 and 3 were 
also elevated in the MSL subtype in comparison to the rest 



of the TNBC subtypes (Additional file 2: Figure SI). Simi- 
larly, analysis of TGFBR3 gene expression across a panel of 
TNBC cell lines, representative of the various subtypes, 
demonstrates that the M and MSL subtypes have relatively 
higher levels of TGFBR3 mRNA (Figure 1C-D). These find- 
ings were validated by qPCR (Figure IE) and immunoblot 
analyses for TpRIII protein levels (Figure IF). Although the 
TNBC M subtype cell lines also showed increased levels of 
TpRIII expression, we focused our studies of this receptor 
on the MSL subtype as their expression is more consistent 
with human datasets (Figure 1A-B). 

Knockdown of TPRIII in MSL TNBC cells leads to 
decreased tumorigenicity in vivo 

In order to determine the significance of the TpRIII expres- 
sion in MSL TNBC cell behavior, we knocked down TpRIII 
in MSL cells and performed orthotopic xenograft tumor 
studies. We used a panel of four shRNA expression vectors 
to optimize TpRIII knockdown, as validated by immunoblot 
and flow cytometry analyses (Figure 2A-C). We utilized im- 
munocompromised nude mice to establish orthotopic xeno- 
graft tumors from cell lines representing the MSL subtype 
of TNBC with and without TpRIII knockdown. Initially we 
tested SUM 159 cells with two shRNA expression vectors 
(TpRIII-KD and TpRIII-KD4) to eliminate off target effects 
of the shRNA (Additional file 2: Figure S2). After establish- 
ing that both expression vectors resulted in a similar pheno- 
type, we used a single shRNA (TpRIII-KD) in all subsequent 
experiments across three MSL cell lines. Knockdown of 
TpRIII in the SUM159 and MDA-MB-231 MSL cell lines 
significantly decreased xenograft tumor growth (Figure 2D- 
E). MDA-MB-157 showed inconsistent results (Additional 
file 2: Figure S3A) and after further investigation we discov- 
ered that the TpRIII-KD tumors expressed TpRIII (Add- 
itional file 2: Figure S3B). Thus, either there was a selection 
against the knockdown in vivo and, therefore, the tumor 
cells expressed TpRIII, or there was a small subpopulation 
of MDA-MB-157 cells at the start of the experiment that 
retained expression and seeded the tumor growth. Regard- 
less, both explanations provide further evidence for the 
tumor-promoting effect of TpRIII. 

Knockdown of TPRIII in MSL cell lines does not affect cell 
proliferation or viability 

Since TpRIII-KD markedly decreased the tumorigenic 
potential of mesenchymal TNBC cells, we further ex- 
plored whether this was due to a proliferation defect. 
TpRIII can bind to all TGF-p ligands but with highest 
affinity for TGF-p2 [29,30]; therefore, cells were treated 
with TGF-p2 in addition to TGF-pl. Both controls and 
TpRIII-KD MSL cell lines responded similarly to the li- 
gands (Figure 3A-B). TpRIII-KD did not alter the prolif- 
eration rates of MSL cell lines (SUM159, MDA-MB-231 
and MDA-MB-157) by live cell counts (Figure 3 A) or 
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Figure 1 TGFBR3 gene expression levels are elevated in the mesenchymal stem-like (MSL) subtype of TNBC. A) Heat map representation 
of gene expression for 587 TNBC tumors for each TGF-(3 ligand and receptor. B) Quantification of average TGFBR3 gene expression across TNBC 
tumor subtypes, average based on individual TNBC tumor probe intensity values (**P = 0.004; ***P <0.0003 for a two-tailed Student's t-test). 
C) Heat map of TGFBR3 mRNA expression in TNBC representative cell lines. D) Quantification of TGFBR3 gene expression across representative 
TNBC cell lines {**P = 0.004; for a two-tailed Student's t-test, ns = not significant). E) qRT-PCR analysis of TGFBR3 average mRNA expression (2" AACt ) 
[28] from representative TNBC cell lines; graph bars represent the mean of three replicates with SEM error bars. F) Immunoblot analysis of T(3RIII 
protein expression in TNBC representative cell lines, results representative of two independent experiments. SEM, standard error of the mean; 
TNBC, triple negative breast cancer; TPRIII, type III transforming growth factor-beta receptor. 



3 H-thymidine incorporation assay (Figure 3B). Consist- 
ent with an intact TGF-p signaling pathway [22,23] we 
have observed an increase in phospho-SMAD2 following 
ligand treatment (Additional file 2: Figure S4). In order 
to examine cell viability and determine whether knock- 
down of TpRIII influenced apoptosis, we analyzed cleaved- 
caspase 3 and cleaved-PARP and we did not detect any 
difference between control and TpRIII-KD MSL cells 
(Figure 3C). 



Knockdown of TPRIII in MSL cells impairs motility, 
invasion and the ability to form invasive protrusions in 
three-dimensional cultures 

Using a validated method (please see methods section 
for details) for measurement of cell migration [21], we 
found that TpRIII-KD significantly decreased the migra- 
tion of SUM159, MDA-MB-231 and MDA-MB-157 cells 
(Figure 4A-C). Treatment with TGF-p ligands had no ef- 
fect on migration. In order to determine the invasive 
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Figure 2 Stable knockdown of TPRIII in MSL TNBC cells reduces xenograft tumor growth. A) Immunoblot analysis of TPRIII protein 
expression in lysates harvested from SUM159 cells stably expressing control and four different TPRIII shRNA vectors (TPRIII-KD, KD2-4). B) Flow 
cytometry analysis of TPRIII protein levels in SUM159 controls, TPRIII-KD and T(3RIII-KD4. C) Immunoblot analysis of TPRIII protein expression in 
lysates harvested from MDA-MB-231 and MDA-MB-157 cells stably expressing control and TPRIII-KD. D-E) Tumors in nude mice were palpable 
three weeks post implantation of the MSL cell line (D, SUM159 and E, MDA-MB-231). Serial tumor volumes (mm 3 ) were measured weekly for five 
weeks post palpation for both controls and TPRIII-KD. Each data bar represents the mean tumor volume of 10 tumors; error bar represents SEM 
(**P <0.005, ***P <0.0001 for a two-tailed Student's t-test). MSL, mesenchymal stem-like; SEM, standard error of the mean; TNBC, triple negative 
breast cancer; TPRIII, type III transforming growth factor-beta receptor. 



properties of MSL lines we analyzed their ability to mi- 
grate through a barrier using an invasion transwell assay. 
TPRIII-KD impaired the ability of the MSL cell lines to 
invade through matrigel pre-coated transwells and the 
addition of TGF-p ligands had little effect on invasion in 
either controls or knockdowns (Figure 4D-F). Next, we 
examined the effect of T(3RIII-KD on the ability of MSL 
cells to form colonies in three-dimensional matrigel 
culture. After five days in culture, SUM 159 cells with 
TpRIII-KD had smooth edges around cell spheres while 
control cells had multiple protrusions invading into the 
surrounding matrix (Figure 4G). These results were 
quantified by calculations of the perimeter, which show 
a significant difference between controls and TpRIII-KD 
(Figure 4H). Overall, these data indicate that TpRIII 
modulates migration and invasion, independent of TGF- 
p stimulation, in MSL cells. To further investigate TGF- 
p pathway signaling [31] in the MSL lines we used 



standard CAGA-luc (Additional file 2: Figure S5A) and 
3TP-lux (Additional file 2: Figure S5C) reporter assays for 
measurement of TGF-p activity [22,23]. Assays were per- 
formed either in the presence of TGF-pl or TGF-p2 ligands 
and compared to untreated controls [30,32]. In addition, 
we performed qPCR analysis for SMAD7 (Additional file 2: 
Figure S5B) and PAI-1 (Additional file 2: Figure S5D) 
gene expression as readout for downstream targets for 
canonical and non-canonical TGF-p activity, respect- 
ively [33,34]. The results of both assays indicate that 
knockdown of TpRIII does not modulate either arm of 
the TGF-p signaling pathway. Thus, MSL lines with 
TpRIII knockdown have resulting phenotypic changes 
without concomitant changes in the TGF-p signaling 
pathways measured. Considering these results and know- 
ing that TpRIII can also bind to bone morphogenetic pro- 
teins (BMPs) [35], we treated the engineered MSL cell 
lines with BMP4. We did not observe significant changes 
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Figure 3 TPRIII-KD in MSL cell lines does not affect cell proliferation. A) Live cell count proliferation assay for SUM159, MDA-MB-231 and 
MDA-MB-157 controls versus TPRIII-KD 72 hours post treatment with TGF-(31 and TGF-(32 ligands; graph bars represent the mean of six replicates 
with SEM error bars (ns = not significant; *P <0.01, ***P <0.0005 for a two-tailed Student's t-test). B) Thymidine incorporation proliferation assay for 
SUM159, MDA-MB-231 and MDA-MB-157 controls versus TPRIII-KD in the presence or absence of TGF-pi and TGF-P2 ligands; graph bars represent 
the mean of six replicates with SEM error bars (***P <0.0004 for a two-tailed Student's t-test). C) Immunoblot analysis of cleaved-caspase 3 and 
PARP protein expression using lysates harvested from SUM159, MDA-MB-231 and MDA-MB-157 stably expressing control and TPRIII-KD. Results 
are representative of three independent experiments. MSL, mesenchymal stem-like; SEM, standard error of the mean;TGF-p, transforming growth 
factor beta; TPRIII-KD, type III transforming growth factor-beta receptor knockdown. 



in Smadl/5/8 phosphorylation in TpRIII-KD versus con- 
trol MSL cells (data not shown). The results suggest that 
TpRIII modulates the tumorigenic potential of MSL 
TNBC cells through other signaling pathways. 



Knockdown of TPRIII is associated with increased 
expression of integrin-a2 in MSL TNBC cells 

To determine which genes and/or signaling pathways are 
significantly altered in MSL cells after TpRIII knockdown, 
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Figure 4 Phenotypic effect of TpRIII knockdown in MSL cells. A-C) Cells were plated around magnetic stencils. After cells had adhered the 
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cell line with TPRIII-KD in comparison to control with or without TGF-(31 and TGF-[32 treatment; graph bars represent the mean of three 
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we performed gene expression microarray analyses on 
SUM159 cells grown in three-dimensional cultures. The 
integrin signaling pathway, along with other cell adhesion 
pathways, were among the most significant pathways 



differentially expressed in TpRIII-KD MSL cells relative to 
control cultures (Additional file 3: Table SI). Analysis of 
individual genes of the integrin pathway revealed that 
ITGA2 was a top gene that was significantly increased 
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upon TpRIII knockdown (Additional file 3: Table S2). 
In vitro qRT-PCR analysis indicates a statistically signifi- 
cant (above two-fold) upregulation of integrin-a2 in the 
TpRIII-KD MSL cells (Figure 5A-C). The upregulation of 
integrin-a2 was further validated by flow analysis across 
all MSL (Figure 5D-F). 

TPRIII modulation of integrin-a2 expression is required 
for the migratory and invasive MSL cell line phenotypes 

Using a clinically relevant, spontaneous mouse model of 
breast cancer progression and metastasis, Ramirez et al. 
demonstrated that integrin-a2pl acts as a tumor suppres- 
sor and a2-null cells were more motile and invasive [36]. 
The in vivo and in vitro findings were further correlated 
with analysis of microarray gene expression datasets of 
human breast and prostate cancers, which showed a 



correlation between decreased expression of ITGA2 and 
poor prognosis. Considering this role of integrin-a2 in 
breast cancer, we hypothesized that the decrease in migra- 
tion and invasion upon TpRIII-KD in MSL cells could be 
rescued by concomitant knockdown of integrin-a2. To 
test our hypothesis, we stably knocked down integrin-a2 
(a2-KD) in the MSL TpRIII-KD cells and performed mi- 
gration and invasion assays (Figure 6A-B and Additional 
file 2: Figure S6A-B). Knockdown of integrin-a2 was suffi- 
cient to reverse the migration (Figure 6C and Additional 
file 2: Figure S6C) and invasion (Figure 6E and Additional 
file 2: Figure S6D) phenotype of MSL cells with TpRIII- 
KD to those of control cells. In addition, using an 
integrin-a2 neutralizing antibody we rescued the migra- 
tory phenotype (Figure 6D) in a manner similar to that 
observed after a2-KD in TpRIII-KD cells. Knelson and 



ITGA2 mRNA Integrin a2 protein 




Figure 5 TPRIII-KD modulates expression of ITGA2 in MSL cells. A-C) qRT-PCR analysis for ITGA2 mRNA expression from MSL TNBC cell lines 
with TPRIII-KD; graph bars represent the mean of three replicates with SEM error bars {**P = 0.003, ***P <0.0007 for a two-tailed Student's t-test). 
D-F) Flow cytometry analysis for integrin-a2 in controls and TPRIII-KD MSL; arrow pointing to the right shows a shift towards an increase in 
protein levels of integrin-a2 in TPRIII-KD MSL cell lines. MSL, mesenchymal stem-like; SEM, standard error of the mean; TPRIII-KD, type III 
transforming growth factor-beta receptor knockdown. 
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Figure 6 Knockdown of integrin-a2 in TPRIII-KD MSL cells reverses migratory and invasive TPRIII-KD phenotypes. A) qRT-PCR analysis for 
ITGA2 mRNA expression before and after stable knockdown of integrin-a2 (a2-KD) in SUM159 cells with TPRIII-KD; graph bars represent the mean 
of the replicates with SEM error bars {***P <0.0001 for a two-tailed Student's t-test). B) Flow cytometry analysis of a2-KD; arrow pointing to the 
left shows a shift towards a decrease in protein amount of integrin-a2 in TPRIII-KD SUM159 cell lines after a2-KD. C) Transwell migration assay 
representing number of cells migrated through transwell upon a2-KD in TPRIII-KD SUM159 cell line (ns = not significant, ***P <0.0001 for a 
two-tailed Student's t-test); bar graph represents a mean of three replicates with SEM error bars. D) Transwell migration assay representing 
number of cells migrated upon treatment of TPRIII-KD SUM159 cell line with anti-a2 blocking antibody (ns = not significant, *P = 0.011, 
***P <0.0001 for a two-tailed Student's t-test). E) Transwell invasion assays with inserts pre-coated with matrigel allowing for testing the number 
of cells that can invade upon a2-KD in TPRIII-KD SUM159 cells (ns = not significant, ***P <0.0001 for a two-tailed Student's t-test). F) Immunoblot 
analysis for phospho-ERK using protein harvested from SUM159 cells with TPRIII-KD and TPRIII-KD /a2-KD. MSL, mesenchymal stem-like; SEM, 
standard error of the mean; TPRIII-KD, type III transforming growth factor-beta receptor knockdown. 



colleagues showed that knockdown of TpRIII leads to 
diminished fibroblast growth factor 2 (FGF2)-mediated 
ERK phosphorylation [37]. Consistent with this previous 
study, after knockdown of TpRIII in the MSL cells, the 
phospho-ERK levels decreased and were restored in the 
cells after simultaneous integrin-a2 and TpRIII knock- 
down (Figure 6F and Additional file 2: Figure S6E). 



Relationship between gene expression of TPRIII and 
integrin-a2 in TNBC patient dataset 

To further investigate the association between TpRIII 
(TGFBR3) and integrin-a2 (ITGA2) in TNBC, we used the 
TNBC patient dataset described in Figure 1A [2] to analyze 
the relationship between TGFBR3 and ITGA2 gene expres- 
sion. Our results indicate an inverse correlation between 
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ITGA2 and TGFBR3 expression across TNBC subtypes. 
In particular, we see the strongest inverse correlation in 
TNBC subtypes with either high TGFBR3 expression 
(MSL; P = 5.274e-06); or low TGFBR3 expression (basal- 
like 2; with P = 5J6e-07 and Luminal AR (LAR); with 
P = 1.759e-07) (Figure 7A-B). The clinical association 
between ITGA2 and TGFBR3 expression is relevant as it 
further links the impact of the interplay between the 
TGF-p and integrins pathways in TNBC. 

Discussion 

Currently, the functional role of TpRIII is controversial 
in breast cancer. Some reports suggest a tumor suppres- 
sive function of TpRIII [18], while other reports indicate 
a tumor-promoting role [19,38-40]. Through GE analysis 
of 587 TNBC patients, we demonstrated that the TGFBR3 
is expressed at a higher level and most consistently in the 
MSL subtype of TNBC. Furthermore, we have identified 
MSL cell line models that express high levels TGFBR3. To 
understand better the molecular basis of TGFBR3 GE 
we used representative MSL cell lines and a TpRIII loss- 



of-function approach. The data presented herein support 
our hypothesis of an oncogenic role for TpRIII in the 
MSL subtype of TNBC. Biologically, knockdown of TpRIII 
in TNBC MSL cell lines resulted in decreased motility 
and invasion, a lack of invasive protrusion in three- 
dimensional culture in vitro and a significant decrease in 
tumor growth in a xenograft mouse model. The observed 
migratory and invasive cell line phenotypes were further 
associated with modulation of the integrin- a2 pathway. 

Previously, the loss of TpRIII expression was corre- 
lated with progression from a pre-invasive to an invasive 
state of breast cancer [18]. In addition, restoring expres- 
sion of TpRIII in a breast cancer cell line led to a de- 
crease of tumor invasiveness in vitro and tumor invasion 
and metastasis in vivo [18]. Other studies have demon- 
strated a frequent loss of TpRIII in human breast can- 
cers [14-17,41]. However, these studies lacked genomic 
analysis of tumors, thus their difference in conclusions 
can be due to the difference in their study models. Tak- 
ing into account the dependency of TGFp signaling on 
the context [42] as well as the heterogeneity of breast 




•o mo mo mo mo mo mo 

m = TGFBR3 O = ITGA2 

Figure 7 TNBC patient dataset shows that expression of TGFBR3 is inversely correlated with expression of ITGA2. A) Heat map 
representation of TGFBR3 and ITGA2 gene expression for 587 TNBC tumors for each TGFBR3 and ITGA2. B) Quantification of gene expression for 
TGFBR3 (solid circles) and ITGA2 (empty circles) across TNBC tumor subtypes using log2 probe intensity values. The whiskers of the boxplot 
extend to the highest value that is within 1.5 interquartile range. P values were generated by performing a paired two-tailed Student's t-test. 
TNBC, triple negative breast cancer. 
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cancers, especially TNBC [2], we took a more focused 
approach to study the role of TpRIII in breast cancer 
progression. Since it has been previously established that 
TpRIII can modulate TGF-p signaling [6,43-46], it is not 
surprising that TpRIII has been shown to have both pro- 
and anti-tumorigenic effects in breast cancer. Our study 
shows that knockdown of TpRIII in MSL cells appears 
not to alter the cells' ability to respond to TGF-p either 
through the canonical and non-canonical arms of the 
pathway, or the BMP pathway. Rather, we observed that 
loss of TpRIII results in a decrease in cell motility and 
invasion. To further investigate potential mechanisms by 
which TpRIII regulates these cell functions, we performed 
gene expression analysis on cells after TpRIII knockdown. 
We found that the expression of genes involved in integrin 
signaling and cell-ECM interactions were significantly dif- 
ferentially regulated after TpRIII knockdown. 

Previous work has shown that inhibition of integrins 
can reverse the transformed state of breast cancer cells 
and that a2pl integrin can play a role in cancer progres- 
sion [47]. A more recent study demonstrated that a2pl in- 
tegrin acts as a metastasis suppressor in breast cancer, 
where migratory and invasive abilities of tumor cells are 
enhanced after loss of a2pl integrin expression [36]. This 
supports our finding wherein a decrease in the migratory 
and invasive phenotype, upon TpRIII knockdown, was 
linked to increased integrin- a2 expression levels. The pre- 
cise mechanistic link between TpRIII and integrin- a2 
expression levels is unknown. The only other association 
between integrins and TpRIII was reported in MCF10A 
breast epithelial cells where TpRIII was shown to regulate 
integrin- a5 localization [48]. 

Knockdown of integrin- a2 in TpRIII-KD MSL TNBC 
cells reverses the loss of motility and invasion that occurs 
upon TpRIII knockdown alone. One explanation for the 
observed rescue of migratory and invasive phenotype is 
through the regulation of ERK phosphorylation possibly 
mediated by integrin- a2. As shown in Figure 6, upon 
knockdown of TpRIII we observed a decrease of phospho- 
ERK simultaneous with an increase in integrin-a2. Fur- 
thermore, upon knockdown of integrin-a2 in TpRIII-KD 
cells we see an increase in phopho-ERK suggesting that 
integrin-a2 is suppressing ERK activity. This is in agree- 
ment with other studies that have shown that integrins 
can regulate ERK activity [49-51]. In addition, studies have 
demonstrated that continuous ERK activity can regulate 
invasion and migration by regulating transcription of 
genes or directly regulating enzymes necessary for cell 
movement [52,53]. Therefore, the increase in phospho- 
ERK seen upon integrin-a2 knockdown could be an ex- 
planation for the increase in mobility of our TpRIII-KD 
cells. Our data show a correlation between TpRIII s modu- 
lation of migration and invasion and the reduction of 
phospho-ERK levels, possibly mediated by integrin-a2. 



Further studies will be required to elucidate the precise 
mechanistic relationship between TpRIII and integrin-a2. 

Conclusions 

In summary, our studies using MSL TNBC models dem- 
onstrate that TpRIII is an oncogenic driver of migration 
and invasion in vitro as well as tumor growth in vivo. Fur- 
ther mechanistic characterization of MSL TNBC would 
provide insights on how to use this protein and/or signal- 
ing pathway as a biomarker or to provide insights to new 
targets for therapy. Considering the limitations of in vitro 
studies, it is necessary to develop a mouse TpRIII breast 
cancer model to further elucidate the role of this mol- 
ecule. Such a model would provide more accurate obser- 
vations for studying the role of TpRIII in the tumor 
microenvironment. The results of this study provide 
mechanistic insight into the role of TpRIII in TNBC and 
highlight an association between TpRIII and integrin-a2 
expression and regulation of cell motility, invasion and 
tumorigenicity. In addition, this study provides a pre- 
clinical rationale for consideration of TpRIII as a potential 
target for further discovery efforts. 

Additional files 



Additional file 1: Additional information about primer sequences, 
antibodies and list of GEO registration numbers [2] are referenced 
in methods section of the manuscript. 

Additional file 2: Figure SI. Average probe intensities for TGF-p 
receptors and ligands across 587 TNBC patients. A-B) Quantification of 
TGFBR1 and TGFBR2 mRNA expression across TNBC tumor subtypes. C-E) 
Quantification of TGFB 1, TGFB2 and TGFB3 mRNA expression. Figure S2. 
Knockdown of TpRIII with two independent shRNA vectors decreases 
orthotopic tumor volume of SUM 159 xenografts. Bars represent mean 
volume of eight tumors. Figure S3. MDA-MB-157 expresses T(3RI 1 1 after 
implanted in vivo thus does not exhibit significant change in tumor 
growth. A) Bars represent mean tumor volume of 10 tumors. B) qRT-PCR 
comparison of TGFBR3 expression in MDA-MB-157 cells before implantation 
and from tumors. Figure S4. pSMAD2 and T(3RII levels indicate that TGF-p 
signaling is intact in TpRIII controls and TPRIII-KD MSL lines. Immunoblot 
analysis. Figure S5. TGF-p signaling appears to remain functional in 
TPRIII-KD MSL cell lines. A) Controls and TPRIII-KD MSL cells were 
examined for CAGA-Luc expression. Bars represents mean of four 
replicates. B) qRT-PCR analysis for SMAD7 mRNA expression; bars 
represent the mean of three replicates. C) 3TP-lux expression. Bars 
represent mean of four replicates. D) qRT-PCR analysis for PAI-1 mRNA 
expression; graph bars represent the mean of three replicates. 
Figure S6. Knockdown of integrin- a2 (a2-KD) in TPRIII-KD MSL cells 
reverses migratory and invasive TpRIII-KD phenotypes. A) qRT-PCR 
analysis. B) Flow cytometry analysis of a2-KD; arrow pointing to the 
left shows a shift towards a decrease in integrin-a2 with TpRIII-KD 
after a2-KD. C) Transwell migration assay with a2-KD in TpRIII-KD; 
bars represents a mean of three replicates. D) Transwell invasion assays 
with inserts pre-coated with matrigel to test for invasion by ct2-KD in 
TPRIII-KD cells. E) Immunoblot analysis for phospho-ERK with TPRIII-KD 
and TpRIII-KD/a2-KD. For all figures, error bars represent SEM, ns = not 
significant and *P = <0.05, **P = <0.01,***P= <0.001. 

Additional file 3: Table SI. Integrin pathway is among significantly 
changed signaling pathways in SUM1 59 TpRIII-KD three-dimensional 
culture system. Genes were considered differentially expressed and 
included for pathway analysis if they met a cutoff of |log2FC| >0.5 and 
FDR adjusted P value <0.05. Pathway analysis was performed by querying 
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against the C2 Canonical Pathways in the Molecular Signature 
Database (MSigDB). Table S2. Integrin family members in SUM159 cells 
three-dimensional cultures with T(3RIII-KD. Table represents list of integrin 
family genes from microarray analysis. Genes are ordered based on 
adjusted P values (low to high). ITGA2 was the top integrin gene with 
lowest P value (P = 0.003). 
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